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^ Abstract 



For a restricted class of potentials (harmonic+Gaussian potentials), we 
express the resolvent integral for the correlation functions of simple traces of 
^ , powers of complex matrices of size N , in term of a determinant; this determi- 

' nant is function of four kernels constructed from the orthogonal polynomials 

corresponding to the potential and from their Cauchy transform. The correla- 
tion functions are a sum of expressions attached to a set of fully packed oriented 
loops configurations; for rotational invariant systems, explicit expressions can 
ly^ ' be written for each configuration and more specifically for the Gaussian poten- 

, tial, we obtain the large N expansion ('t Hooft expansion) and the so-called 

BMN limit. 



1 Introduction 



^ I Given a potential V {M, M+) where M is a x matrix with complex 

elements, we define the partition function as 

Z = JdM dM+ e-^'- (1.1) 

where we integrate over all matrix elements. We define the correlation functions 
as 

< []0, (M, M+) >- ^ / dM dM+ []0, (M, M+) e"^"^ v{m,m^) (1,2) 



Most observables O (M, M+) can be expressed in terms of traces of product of 
the matrices M and M"*"; among these traces, we distinguish between "mixed 



traces" like Tr M-'^ (M+)^' M-'^ {M+ f^ 



which are also called " words" , 
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and simple traces where there is only one type of matrix {M or M'^) inside the 
trace. In this publication, we restrict ourselves to the case of simple traces 



p q 

< JJ [Tr M-^'] [| \ Tr (M+)^*] > and Ki>Q (1.3) 

i=l i=l 

for a restricted class of potentials (described below). We now explain the method 
used to obtain the above correlation functions. 

If certain conditions on the potential V (Af, Af+) arc satisfied, we were able 
in Ref.[l] to calculate the following expectation values 

Hdet (r;*-M+) ndet(^.-M) 

\ I dM dM+ '-^ l^i e-^^ v{MM+) (1.4) 

JJdet«-M+) Y[Aet{yi-M) 

i=l i=l 

for any numbers Li, Mi, L2, M2 of determinants. The general result is essen- 
tially a determinant constructed from four kernels, each kernel is a function of 
the external sources ^i, ry*, yi, x* which enter the orthogonal polynomials rela- 
tive to the potential V (z, z*) and their Cauchy transform. The structure of this 
determinant is very dependant of the relative values of Li, Mi, L2, In this 
publication we are only interested in the special case where Li = Mi, L2 = M2 

M2 Ml 

since in that case, if we perform the derivatives (off) (^) 
(1.4) we obtain 



, ^ M2 Ml 

\ / dM dM+ r]Tr—^ TTrr— ^ e-^^ v{m,m+) (1.5) 



The above resolvent integrals are the generating functionals for the integrals 

. M2 Ml 

^ dM dM+ Y[Tr [(M+)'^-] []Tr [(M)^*] e"^'' v(m,m+) ^^ ^-^ 

•' i=l i=l 

which describe the correlation functions (1.3). 

Let us now describe the conditions required by the potential y (M, M+) . 
Any matrix M can be triangularizcd by a unitary matrix 

M=U [D + T] U+ (1.7) 

where D is the diagonal matrix of the eigenvalues {zi,...,zn) of M, and T is an 
upper strictly triangular matrix. Consequently, any power of M and M+ is of 



2 



the form 



M'' = U [D'' + Tk] U+ 



U+ 



(1.8.a) 
(1.8.b) 



where Tk is an upper triangular matrix which also depends of the eigenvalues 
for fc > 1. Clearly, we have 



N 



i=l 



i=l 

N 

Tr {M+f = Y,{z*f 

On the other hand, 

Tr MM+ = ^ z* + Tr TT+ 

i=l 

In the transformation (1.7), 

det (Ci -M) = det (Ci - D) 
det {ri* - M+) = det (r?* - D+) 

and the integration measure becomes 

dM dM+ = dU dT dT+ dD dD+ ]J 



(1.9.a) 
(1.9.b) 



JV 



(1.10) 



Zi — Zi 



(l.ll.a) 
(l.ll.b) 



(1.12) 



We note that dM dM^ represents 2N'^ variables of integration while dU represents 

(TV — 1) variables of integration, dT dT~^ represents TV {N — 1) variables of 
integration and dD dD^ represents 27V variables of integration. The Jacobian 
of the transformation is given by the product of two Vandermonde determinants 
A{z) A{z*) where 



A{z) = l[{zi-Zj) = 



T^N-l (zi) ... TTJV-I (zn) 

7ro(2i) ... 7ro(zjv) 



(1.13) 



and where the polynomials {ttq (z) , tti (z) , 7r2 {z) , ttjv-i {z) }are any set of 
monic polynomials of successive degree 0,1,2,...,A'' — 1. 

Consequently, if wc choose the potential V (M, M+) to be of the form 

V {M,M+) = MM+ + V{M) + V{M+) (1.14) 
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wc may separate the integrations over the unitary group dU and over the upper 
triangular part dT dT~^ e~'^^ and these two contributions cancel in the 
ratio of the numerator and the denominator of (1.6). We are left with integrals 
over the eigenvalues only. 

More generally, we consider the integrals Zff, Ipf and Jjv and a potential 
V {z, z*) such that the following integrals exist 



= j dn{z,z*) 



(l.lS.a) 



In = 



1 r ^ 


r i2 

n - 

i=l 


-V*) n(^i 

i=l 


M2 
Li=i 


Ml 
i=l 




1 Ml ( N 

.HE 


N 

Y[dzi dz* l^i - Zj 

i=l i<j 


e 





Vi) 



The integrals Jjv are generating functionals for the integrals 



^ . M2 ( N \ Ml ( N 

^fd,iz,zn n E(^r n E 



(l.lS.b) 

1 

(1.15.C) 
(1.15.d) 

(1.16) 



A second condition on V {M, M+) is that V {z, z*) is real and admits a set 
of orthogonal polynomials so that the technique of Ref[l] can be applied (the 
reality of V (z, z*) could be forgotten at the price of introducing biorthogonal 
polynomials but we shall not consider this case). We introduce the infinite set 
{Pn {z)} of orthogonal, monic polynomials such that 



Pi {Z) Pn {Z) e-^(^'^*) = K Snm (1.17) 

where p^ (z) is a short notation for [p^ {z)]* ■ Then, we consider the measure 

i[{z*-v*) n(^-^*) 

i=l i=l 

M2 Ml 



dn{z,z*;ii,7j*;y„x*)=d'z'-^ e-^(^'^*) (1.18) 
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Wc proved in Rcf.[l] the existence and we constructed explieitely the set of 
monic, biorthogonal polynomials g„ {z;^i,ri^;yi,x1) and g* {z;r]i,S,l;Xi,yl) sat- 
isfying 

j dn{z,z*;^i,r]*;yi,x*) q^{z;r]i,^*;Xi,y*) Qn iz;Ci,r]*;yi,x*) = \\qnf 5nm 

(1.19) 

This construction is a generalization of Christoffel's result [2] which shows 
that given a positive Borel measure of one variable fi// (x) on the real line and 
its infinite set of orthogonal polynomials, it is possible to construct an infinite 
set of orthogonal polynomials for the measure 

L 

dii{x;^i)=dii{x) Jl{x-^i) (1.20) 

i=l 

This result was extended to measures with external sources at the denominator 
by Uvarov [3] in 1969 and recently by Fyodorov and Strahov [4] . In 2003, 
Akemann and Vernizzi extended ChristofFers result to measures on the complex 
plane with external sources at the numerator [5] (see also Ref.[6]); their work 
was generalized to sources at the denominator as well in Ref.[l] and [7]. 

The consequence of the existence of biorthogonal polynomials for the mea- 
sure (1.18) is that 

JV-l 

Zn = N\ ]J hi (1.21.a) 

i=0 

In = N\ [] llftf (1.21.b) 

In the case L\ = Mi, L2 = M2, the pseudonorms \\qi \f are found in Ref.[l] to 
be the ratio of two determinants 

hif = h^ i > (1.22) 

J-'i-l 

where -D„ is the determinant 



Nn{^i,yi) 

Nn (Cl,2/Mi) 



Nn{iMr,yi) 
Nn {S,Mi , yMi ) 



An {xt,yi) 
N* {VM2,xi) 



An {x*M^,yMi) 

K{Vi,xm2) 

N* {VM2,XM2) 

(1.23) 
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and _D_i is defined in (1.32). Consequently, in the case Li = Mi, L2 = M2 
the integral In is found to be 



Mi(Mi-l) MoCMo-l) 

In = (-)^^ 



M2 

n - '^t, 



Ml 

n (y* - 



A(a;*) A(y) A (77*) A 



D 



N-l 



(1.24) 



M2 



Ml 



Then, by apphcation of the derivatives ( j n(^) on (1.24), 

we proved in Ref.[l] that 

Jn = "Dn-i'' (1.25) 
where is the "subtracted determinant" 

ifn(a,Cl) - A^n(CMi,ei) ^«(??^,Cl) - An{v*M,,^l) 

^n{^l,S,Mi) ■■■ Hn{£.Mn£,Mi) (j?! : ^Mi ) ••• ('7m2 ' '^^^i ) 

Kn{ii,vl) ■■■ Kn{iMi,vl) Hnim^m) ■■■ N*{'qM2,m) 

^n(Cl,??M2) - Kn{^M,,V*M2) ^ninum^) - ^^n (»?M2 , »?M2 ) 

(1.26) 

The determinant Z)„ is expressed in terms of four kernels Kn, Nn,N*, An 
which are defined with some of their properties at the end of this introduction. 
The determinant is obtained from the determinant Dn by changing all 

kernels Nn on the diagonal into i?„, changing all yk into S^k, all x*/, into i]l and 
finally by ignoring all double poles at = ^fc and at rj* = -ql (this opera- 
tion is denoted by " ") as we develop the determinant (1.26) with terms like 
Nn {^j,^k) (a,$j) or TV,: {r)k,r]j) N* {Vj,Vk)- In fact, from (1.30.b) we ob- 
serve that there is no single poles either at £,j = ^k or at ij* = r/l since the 
residues are zero 



'iv„(c„a) Tv„(a,e, 



Hn (Cj ) Cfc ) Hn {^k , Q ) 



^3 - Cfe 



(1.27) 



In section 2, we develop the determinant in terms of fully packed 

oriented loops configurations. Then, in section 3, we give an expression for each 
configuration in the case of rotationnally invariant systems ; finally in section 4, 
we use the Gaussian potential to obtain the correlation functions (1.3) which in 
that case are nothing but the number of graphs obtained from (1.6) by Wick's 
contraction. The large A'' expansion (where N is the size of the matrix) is 
considered and provides the number of graphs according to their genus (t'Hooft 
expansion [8]); another limit is also considered, namely the so-called BMN limit 
[9], where J^, Ki in (1.3) and A'' are large but the ratios and are constant. 
Similar results were presented in [10] and [11] . 
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The end of this introduction is devoted to the description of the kernels 
Kn, Nn, An and to some of their properties. We refer to Ref.[l] for the technical 
details; here, we simply give the definition of these kernels 



i=0 



P» (0 P*i iv) 
hi 



(1.28) 



We introduce the following integrals (in the case of potentials of one variable, 
the integral (1.29. a) is called the Cauchy-Hilbert transform of the corresponding 
polynomial) 



then. 



tn {v) = ( d'^z pI {z) e" 

J z-y 

Q{x*,y)= I (fz -— ^ r 

J {z*-x*){z-y) 



1=0 



V{z,z*) 



-V{z,z*) 



1 



Hni^,y) 



i=0 



hi 



Q{x*,y) 



In (1.22), we also use the extensions 

N-i{^,y) = 



y-^ 

A-i{x*,y) = -Q{x*,y) 



and the determinant 



A/-_i(0,?a) A^iix*,y,) 
Nl^{7]k,xi) 



D-i = 

Z)_i = |iV_i(^„yi)l • \Nl^{ilk,xi) 



where 



1 



1 



|iV_i(e„2/i)| = (-) 



A(t/)A(0 



Jl{yi-^j) 

hi 



(1.29.a 
(1.29.b 

(1.30.a 
(1.30.b 
(1.30.C 

(l.Sl.a 
(l.Sl.b 
(1.31.C 

(1.32.a 
(1.32.b 

(1.33.a 
(1.33.b 
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In Ref.[l], we prove the following properties 



E 

i=0 



hi ^-y 



(1.34.a) 



i=0 

so that we have the formal power series 

PiiO U{y) 



i=n+l 

Finally, we have the asymptotic behaviour 



An{x*,y) = - £ telAM (i.35.b) 



tn{y) asy^oo (1.36.a) 

^"(^'2^)^^^^^^ asy^oo (1.36.b) 
^^^^ asy^oo (1.36.C) 



2 The graph description of the determinant " 



The determinant "-D„" as given in (1.26) is a sum of products of oriented prop- 
agators joining the various points (for simplicity, we label the points and the 

corresponding variables by the same letter). Let us split the plane into two parts 
separated by a border line; on the right part we draw the various points rj* and 
in the left part we draw the various points ^j. As we develop the determinant, 
every time we meet a kernel An (ViT^.j) we draw an oriented propagator from 
the point 77* towards the point crossing the border from right to left; similarly 
with the kernel Kn {^1,11 j) we draw an oriented propagator from the point 
towards the point rj* crossing the border from left to right. Then, inside the 
^-region we draw an oriented propagator from the point towards the point 

corresponding to the kernel Nn {^i, in addition, the diagonal part of 
defines self closing loops on one point from the kernel _ff„ (^i,Ci)- Similarly, 
inside the 77*-region we draw an oriented propagator from the point rj* towards 
the point rj* corresponding to the kernel N* (jij , rji ) : in addition, the diagonal 
part of "£)„" defines self closing loops on one point rj^ from the kernel H* [rji, jji). 
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As wc develop the determinant, the product of propagators describes a num- 
ber of disjoint oriented closed loops 7a (including self-loop with one point) so 
that any point or r]* belongs to one closed loop 7a and only one; we call such 
a configuration " fully packed oriented loops configuration" which we denote by 
zu (fig.l); the number of oriented loops in zu is denoted # {zu) 

^ = ^tifja (2.1) 




fig.l : contribution to Wg.e {^i, . . . ,^s;Vi, ■ ■ ■ , Ve) 



The determinant "-D„" is the sum over the (Mi + M2)! configurations w 

"Z)„" = (-)^^+^^ J2 n\-/.J (2.2) 

a=l 

Given an oriented loop 7, which is not made of N propagators alone and/or of 
A^* propagators alone, the amplitude /-y is the product of a certain number g > 
of propagators Kn and of the same number q of propagators A„; the remaining 
propagators describe q chains of propagators A'^ and q chains of propagators 
A^* (such chains may be reduced to a single point). We call skeleton graph an 
oriented loop 7 without N and N* propagators. To any oriented loop 7 we 
associate; a skeleton graph by contraction of each of its chains of propagators A^ 
or A"* into a point. 

Along a given chain, the propagators A' or A'* have poles when the ^'s or the 
ri's coincide. However, these poles cancel when we sum over all permutations of 
the ^'s or of the rj's inside a given chain of propagators A^ or A^*. In appendix 
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A wc show how the poles cancel when we sum over several configurations which 
have the same associated skeleton graph. 

For a given potential V {z,z*) , we must calculate the orthogonal polyno- 
mials Pn (z) and their Cauchy transform t„ (z) in order to evaluate the kernels 
which are the propagators of the graphs. Of course, each case is a specific case; 
however, some general structure can be obtained in the special case where the 
system (potential and domain of integration) is rotationnally invariant. 



3 The |, ^ expansion in the case of the rota- 
tionnally invariant systems 

The expressions for the four kernels Kn , Nn , N* , An simplify greatly in the case 
of rotationnally invariant systems because the angular integrations in (1.17) is 
trivially performed. Given a potential V (zz*) and a disc of radius R which 
might be the complex plane (i? — oo) if cxp {—V (zz*)) decreases strongly 
enough at oo. In that case, the orthogonal polynomials are the so-called Gini- 
bre's polynomials 

Pn {Z) = (3.1) 

and their norms squared are 

hn = 2TT / dp e-^ip') (3.2) 

Jo 

In the Gaussian case where V {zz*) = zz* and R = oo, 

hn=nr{n+l) (3.3) 
The function t„ (y) defined in (1.29. a) is found to be 

if R = oo this result is still valid up to exponentially small terms when 
and can be used as a formal power series (see also (1.34. a)). 



According to (1.28) and (1.30) the formal power series for the kernels 
Kn, Hn,Nn and An are 

Kn{^,V*) = ^^-^ (3.5.a) 



, rrn+l _ tn+l 



Hn {^uii) = --^ ' , J ^ (3.5.b) 
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so that 



77-1-1 

Hn{^i,^i) = ^ (3.6.a) 



Finally, 



A^nte,^,)= (gj ^-3^ (3.6.b) 



7 

A„(^*,0 = - E 77T7TT (3-7) 



It is now possible to calculate for any closed loop 7 the corresponding am- 
plitude I J which contributes to the value of the determinant "Dn" . 



1°) the |, expansion for the skeleton graphs: given a skeleton graph 
with g > propagators A and q propagators K, the expression for 7^ can be 
written 

Ij = f[Ar,{v*^^,U) flKn{U,Va,+,) , with = 7?*, (3.8) 

fe=l i 

As we expand the propagators A„ and JC„ according to (3. 5. a) and (3.7) we 
introduce an internal momentum for An (r?^^. , S,ak ) ^'Hd. an internal momentum 

jk for K„ ^^ctfc , Vak+ij ' summation over ik runs from n + 1 to oo while the 
summation over jk runs from to n. By convention, we extend the definition 
of hj to negative j's as /ij<o = oo so that in Kn we may eventually ignore the 
lower limit of summation. 

The expansion of 7-y is of the form 

A= E 11(7^1) n(;i^)* nJ^,Kk) (3.9) 

{J's>0,K's>0} k=l ^?"'= ^ fe=l ^ 

where 

/(J.,i^.) = (-)^ E (3-10) 

and D{jk.Kk} is the summation domain for the indices ife and jfe- The determi- 
nation of this domain is the main difficulty. Clearly, all the indices i' s and j' s 
are fixed in terms of the J's and the K' s up to a translation. 

We first describe some simple cxemples: 
a) let us consider the graph 71 (fig. 2). 
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fig. 2 : graph 71 W^^ {x, y) 
In that case, the momentum conservation law is 

i-j = J = K (3.11) 

so that the domain of summation for the momentum j is 

n- J + 1 < j <n (3.12) 

Consequently, 

n , 

j=n-J+l ^ 

The equation (3.13) is still valid in the case where n — J + 1 < with the 
convention /ij<o = 00 ■ We define the characteristic function for the skeleton 
graph 7i 

n , 

$,,(n,J)= (3-14) 

then, 

I^, ( J, K) = -5j,K (1 - Tj) (n, J) (3.15) 
where the operator Tj is a translation operator over n 

Tj f{n) = f{n-J) (3.16) 



12 



b) let us consider the graph 72 (fig. 3). 



fig. 3 : graph 72 — > contribution to ^2,2 (2^1, 2:2; j/i, 2/2) 

Choosing j2 = j as the independant loop momentum, we observe the follow- 
ing constraints 



j >n + l- Ki, j>n+l-J2 

j < n — Ki + Ji, j < n 
J = Ji + J2 = Ki+K2 = K 



(3.17.a) 

(3.17.b) 
(3.17.C) 



Consequently, we obtain 



Ij2 {Ji, J2, Ki, K2) = Sj^K O {Ji, Ki) (n, Ji, Ki) 

O {Ji, Ki) = Tsup{Ki-Ji,0) - Tinf{Ki,J2) 



^72 {n,Ji,Ki) = 



— ' hj+Ki-ji hj 



(3.18.a) 
(3.18.b) 

(3.18.C) 



c) the generalization to the skeleton graph 7, with q pairs of propagators An 
and Kn is now straightforward. We define 

J = ^Jfe K = J2Kk (3.19) 
fe=i fe=i 

and we choose j = jg as independant loop momentum so that 

ik = jq + bk jk = jq + ak (3.20) 

with 

k 

ttk = 'Y{Kp- Jp) , a,=o (3.21.a) 

p=i 

bk=ak + Jk>l (3.21.b) 
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Wc note from (3.20) that any variable bk is strictly larger than any variable ak'- 
The contribution of the graph 7, is 

J^, {Ji, Ki) = {-Y 6j,K (Tx - Ty) (n, Ji, K^) (3.22.a) 
(n, J„ K,) = Y: '^'^ (3.22.b) 

j=0 '^i+'^i ■■■■'^j+aq 



with 



< X = Sup (flfc) (3.23.a) 
1<Y = Inf {bk) (3.23.b) 



2°) the |, i expansion for any graph: a one loop graph 7 is made 
of q propagators A, q propagators K, q drains (Ci,...,Cq) of propagators N 
and q chains (CJ',...,C*) of propagators N* (the chains can be reduced to a 
single point). The case q = in made of loops of propagators alone, or 
of N* propagators alone. To each chain C'i we associate an external outgoing 
momentum 

Ai= J2 (3.24) 
and to each chain C* an external ingoing momentum 

e. = ^ Kk (3.25) 

Now, to each graph 7 we associate a skeleton graph s (7) obtained from 7 by 

shrinking into a point each of the chain Ci or C* . To each of these points we 
associate the corresponding external momentum Aj or O^; we have 

Q q 

J = Y,\ = Y.^i=Y.Jk=Y.Kk = K (3.26) 

Because all J's and K's are positive or nul, equation (3.26) shows that if the 
graph 7 is made of propagators alone {q = 0), since the K's are zero, the J's 
are also zero. The only non zero contribution in this category is the self-closed 
loop with one point which is Hn {^i, ^i) = — The two points loop is 

"Ni^u^j) N{^j,^i)" =0 (3.27) 

Larger loops of propagators alone, or of N* propagators alone are shown in 
appendix A to be zero. 

The functions hi are only present in the A and K propagators so that the 
function to be summed for a graph 7 is the characteristic function of the skeleton 
graph s (7) 

^^•+''^-^^'+'"' (3.28) 
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where j is the momentum corresponding to the propagator Kg and 

k 

p=i 

&fe = flfc + Afe > 1 



(3.29.a) 
(3.29.b) 



(3.30) 



The equation (3.10) for a one loop graph 7 with propagators A'' and A''* can be 
written 

j-Q "■j+ai-'-il'j+aq 

where the operator 0{Ji,Ki) determins the range of summation for the variable 

j and is described below. 

We introduce the commutative operation (8) 



Ta®Th= Tsupia.b) 

1 (8) Tb = Tsup(i[),b) 



(3.31.a) 
(3.31.b) 



so that conditions like n — a + 1 < j < n together with n — 6 + 1 < j < n, where 
a and b are positive, can be written 



J2 ... = (i-T„)®(i-r,) 

j=n-Inf{a,b) + l j=0 



(3.32) 



Now, we proved in appendix A that for a given chain d with p points ((p — 1) 
propagators), the sum over the p\ permutations generates a domain of summa- 
tion for the variable ji described by the operator JJ^ (1 — Tj^) ; since the prop- 

agator with momentum ji is between the chains Cj and C*_^_i (with Cg+i = C*) 
the constraints on ji can be written 



n (^-tkj 



(3.33) 



and using the relation (3.20) we can write the corresponding constraints on jg as 



T 

-'■ ai 



n (1-^^^) 



n (^-tkj 



(3.34) 



We collect all the constraints on the various momenta ji for i = l,...,q and 
write 



o(Ji,/^0 = ®Li 



Ta, 



n (^-Tj.) 



n (^-tkj 



(3.35) 
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4 The Gaussian potential, large N and BMN 



expansion 



We now apply the results of section 3 to the Gaussian potential with the function 
hn given in (3.3); we also replace the index n of section 3 by {N — 1) according 
to (1.25). 

1°) let us consider the two points function 



Wi,i iN,x,y) = <Tr- 



1 



Tr- 



1 



X - M y- M+ 
which is equal from (1.26) to the determinant 



> 



N-l 



det" 



Hn-i{x,x) AN-iiy,x) 
KN-i{x,y) HN-i{y,y) 



(4.1) 



(4.2) 



The part Hjq^i {x,x) H^^i {y,y) of the determinant corresponds to a discon- 
nected graph made of two self closing loops and from (3. 6. a) is equal to 
The part Ajv-i {y, x) K^-i {x, y) of the determinant corresponds to the graph 
7i (fig. 2). The corresponding function {N — 1, J) is given in (3.14) 



JV-l 



(7V-1,J)=E^ 



r (j + J + 1) 1 r (iv + J + 1) 



i=o 



+ 



J + i 



r(Ar) 



(4.3) 



Consequently, according to (3.15) 



An-\ {y, x) Kn-i {x, y) = - ^ 



\{xyy+' J+l 



r(7v + j + i) r(iv + i) 



T{N) 



The square bracket [ ] in (4.4) has clearly a expansion 



T{N - J) 
(4.4) 



r(Ar + J+i) r(7V + i) 



T{N) 



T{N-J) 



= 2 N' 



J 



E 

p odd>0 



NP- 



(4.5) 



where ap (J) is described in appendix B. After a convenient rescaling, we write 
AN-i(yVN,xVN^ KN-i(xVN,yVN^ as 



1 °° 1 



J + 



24iV2 



+ (j_2) + 0(^ 



1 

"iV 



1 xy {3xy + 2) 



{xy-lf iV2 {xy-lf 



(4.6.a) 
(4.6.b) 
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The so-called BMN limit of J^, {J,K) (defined in (3.15), (4.3)) is the large 
N, large J limit with constant j^, ^ ■ From the large J behaviour of Up (J) 
as given in appendix B, we obtain the BMN limit as 

I^AJ,K)^-6j,K — «/^(^) (4-7) 



At this stage, we find convenient to introduce a formalism developped in 
appendix B and which makes the expansion and the BMN limit transparant. 
We define 



J 



J 



^«(j+i) ^ \ JV-- + 1... A/'+- 



J 



(4.8) 



which has manifestly a -^rr^ expansion 



^"(J+i) = J2 {-)" Sap {J) nJ-^p+^ 

p=0 



(4.9) 



where the coefficients T,2p (J) are given in appendix B. Consequently, the ex- 
pression (4.4) can be written as 



1 T{N + J+1) 



J+1 



T{N) 



7^1 ^) ^'-''^ 



which is manifestly a expansion. Also the BMN fimit of iV**(-^+i) is simply 
N-^+^ so that the BMN fimit of (4.10) consists in replacing ^ by and 

_/V**(-^+i) by iV''+^ 



2°) we now consider the three points function 



W2,i {N,xi,X2,y) = <Tr- 



1 



Tr- 



1 



1 



xi-M X2- M y- M+ 
which is equal from (1.26) to the determinant 



> 



JV-l 



det" 



Hn-i{xi,xi) Nn-i{x2,xi) AN-i{y,xi) 

Nm-i{xi,X2) Hm-i(x2,X2) AN^l(y,X2) 

KN-i{xi,y) KN-i{x2,y) HN-i{y,y) 



(4.11) 



(4.12) 



The connected contribution to (4.12) is 

[AN-i{y,xi) Nn-i{xi,X2) Kn-i {x2,y) + {xi X2)] ■ 
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fig. 4 : contribution to WJ^ {xi,X2;y) 
In section 3 and in appendix A, we prove that this combination is equal to 



111 

^2,1 W y)= ^TT ^2 -^1' "^2, K) 

X-\ ^ 



/2 (Af, Ji, J2, AT) = (1 - TjJ (1 - TjJ 



(4.13.a) 



1 r(7V + J+l) 



j + 1 r(iv) 



(4.13.b) 



(4.13.C) 



where J = K = Ji + J2. The large TV behaviour of 
W|,i (A,a;iViV,X2//V,y//v) is found to be 



y s—^ 1 1 

^+24^^"^"^^ (J-2){(J+1) {J-2)-2JiJ2} + o(^ 



or 



Wa'i (N,xis/N,X2\fN,yy/N^ = - 



2y_ 



y^xiX2 - 1 



+ 



_{yxi - if {yx2 - if 
Finally from (4.13.c), the BMN limit, where N, Ji, J2 are large with 
-j^ constant, is obtained as 

l2{N,J^,J2,K)r^-6j,K ^^[w 



(4.14) 



1 

(4.15) 



(4.16) 
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3°) the generalisation of the situation 2°) is the calcumlation of 
Wp i {N,xi...,Xp,y) and is easy to understand. We sum over j?! graphs ob- 
tained by exchanging the p points Xi and we calculate 



w;^,{N,Xi,y) 



p-i 



AN_i{y,xi) Y\_^N-i{xi,Xi+i) KN_i{xp,y) + Sym{x., 

1=1 

(4.17) 




fig. 5 : contribution to Wpi {xi,. . . , Xp-, y) 

From section 3 and appendix A, we know that the poles of the propagators 
Nn-1 {xi, Xi-i-i) cancel in the sum and that the result can be written as 

oo p ^ ^ 

{N, Xi, y)= J2 U 7j<TT ^) (^-l^-^) 

{Ji, if}=li=l ^ 



Ip{N,J,,K) = -5j,K l[{l-Tj^) -j—^ 



1 T{N + J+1) 



T{N) 



2P 



Ip {N, Ji,K) = -Sj,K fj^U'^i^^^ 



i=l 



J, _d_ 

2 dN 



(4.18.b) 
(4.18.C) 



where J = K = Ji + ... + Jp. The asymptotic behaviour at large N of 
W (^N,XiVN,yVN^ , after a proper rescaling of the variables, is 



p_l oo p 



Ji+1 



Ji...JpJ{J -1)...{J -p + 2)+0 



(4.19) 
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or 



3p-l 



N 



2p-2 



i=l 



o 



Ar2 



The BMN limit is trivially obtained as 



IpiN,J,,K) ~ -6j,K 



1 P 



1=1 ^ 



(4.20) 

(4.21) 



4°) We now apply the formalism to the four points, connected, resolvant 
function W2 2 {N, xi,X2,yi,y2)- First we calculate the contribution of the graph 
of fig. 6 and its crossed symmetric ones (Ji J2 and Ki K2). 



Jl ^ XI 



fig. 6 : contribution to W22 {xi,X2;yi,y2) 



We write 



W{N,x^,x2,yi,y2)= £ (n -JT+t) (n 

{Ji, K}=1 \i=l ) \i=\ J 

(4.22) 

The formalism of section 3 and appendix A applied to the corresponding four 
graphs gives 



1 r(iv + ./ + i) 



J+l T{N) 

0{J,,Ki) = [(1 - T,,) (1 - Tj,) ® (1 - Tk,) (1 - Tk,)\ (4.23.b) 
0{Ji,Ki) = {l-Ti^f){l + Tsup) 

where the symbol (8) is defined in (3.31) and where 



(4.23. a) 

(4.23.b) 
(4.23.C) 



J ^ K = Jl + J2 = Ki + K2 (4.24.a) 
Inf = InfiJi,Ki) Sup = Sup{Ji,Ki) (4.24.b) 
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The function I {N, Ji,Ki) can also be written 



I(N,Ji,Ki) = -SjK^- sh ( LulA. 

The BMN Umit is immediately obtained as 

r/^r . . 4Af''+i , fJ.Inf\ , fJ.Sup 

I (N, J„ K,) ~ -Sj,K^ sh j ch 

The largo iV limit is obtained by expansion of the sh and ch in (4.25) as 
I{N,J„K,) = -dj,K2Inf N'^ * 



(4.26) 



(4.27) 



5°) The remaining part of the connected four points resolvant function 
W22 (^; xi,X2,yi, 2/2) is attached to the skeleton graphs of fig. 3 and its crossed 
symmetric one ( Ji J2) or {Ki K2) . According to section 3, we have for 
the graph of fig. 3 



JV-l 



h {N, Ji, Ki) = Sj,K {Tx - Ty) ^ 



r(j + i^i + i) r(.? + J2 + i) 



^ r(j + i^i- Ji + i) T{j + i) 

with J = + J2 = Ki + K2 = K and with 

X = Sup{Kx- Ji,Q) Y = Inf{Kx,J2) 

The sum in (4.28) is transformed in appendix B (B.31) into 

r(iv + j-p + 1) 



Inf(Ji,K2) 
p-0 



T{N) 



(4.28) 
(4.29) 

(4.30) 



We can now compute the large behaviour of /i {N, Ji,Ki); for instance, in 
the sector Ji >K\,K2 > J2 we obtain 



h{N,J,,K,) = 6j^K J2 



(4.31.a) 

r = J (J - 1) [3 + J| - ( J + 2)] - 12 K1K2 Ji ( Ji - 1) 

(4.31.b) 

For the BMN behaviour, each term in (4.30) contributes, however the sum of 
these contributions can be performed and we obtain 

h (iV, JuKi) ^ 5j,K — ^ '^\^] (4.32) 
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If wc sum the graph of fig. 3 and the crossed one, and if we describe each of 
the sectors as 

Sup >A,B> Inf (4.33) 

with A + B = Sup + Inf = J = K, we obtain the large N behaviour for both 
graphs as 



I{N,Ji,Ki)=Sj,K2Inf 



(4.34.a) 



r = J ( J - 1) [3 Sup^ + Inf^ - ( J + 2)] - 12 AB Sup {Sup - 1) 

(4.34.b) 

and the BMN behaviour as 

nM,J,.K,-,^>,. '-^ ,, (^),.((^^^) (4,35) 

6°) Finally, we sum the contributions 4°) and 5°) to obtain the asymptotic 
behaviour for the complete connected four points correlation function; taking 
into account the relative minus sign between 4°) and 5°), we obtain 

r§ (AT, Ji, Ki) = -Sj,K N'-^ J1J2K1K2 {Sup - 1) + O {N''-^) 

(4.36.a) 

rc(M T A 8iVJ+i ,fJ.Inf\ fA.Sup\ [B.Sup 
{N, Ji, Ki) ~ -6j^K J — sh — — sh sh 



J \ 2N J \ 2N J \ 2N 

(4.36.b) 



We close this section with the following remarks. Although, in principle, we 
should be able to calculate the contribution of any graph, wc arc unable at the 
present time to give the expression for the sum of the graphs which enters a given 
fc-points correlation function; this situation may improve if we use the so-called 
loop equations[12] . From the exemple 5°), we sec that the t'Hooft topological 
expansion in is not true graph by graph and that the cancellations between 
graphs are far from evident. Finally, we verify on these exemples that at large 
N the resolvant function properly rescaled behave as 

Wk„k, (N,Xiy/N,y^VN) ^ O [n^) (4.37) 

because of the disconnected one point self closing loops (the lowest term of the 
expansion) and that the connected resolvant function properly rescaled behave 
as 

W^^^^^ (iV, Xi^/N, Vi^/N) ^ O (j^) (4-38) 



where k = ki + k2- 
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A Appendix: summation over the chains of prop- 
agators N and A^* 

We consider a chain of {p — I) propagators N (the arguments are the same 
for the chains of propagators N*). The propagator (6, 6+1) has a pole at 

= 6+1. However, when wo sum over the p\ chains obtained by exchanging 
the p points, the poles disappear as we explain now: we consider two points 

and and we look for the pole at ^5 = ; many chains do not contain 
the propagators iV„(^h,^c) or Nn{£,c,^b) and have no pole at = ^c- Now, 
we associate by pair the chains with a propagator iV„ (^6,^c) or N„ {^c,^b) by 
exchanging b and c; the combination 

{Nn{^a,^b) Nni^b,^:) Na{^,,^d)+Nn{^a,^c) A^n(6,6) (6, ^d)) haS UO 

pole at £,b = £,c and is equal to 



Similar arguments prove the absence of poles at ^b = if ^ or/and c are end- 
points of the chain, for instance in{Nn (^6, ^c) + (Cci ^b)) or in 
(Nn (Ca, 6) N„ i^h, ^c) + N„ (Ca, Cc) (Cc, 6)) • We remind that the self clos- 
ing loops Nn{^bj^b) are replaced in (1.26) by Hni^b^^b) given in (3. 6. a) and 
contribute only to disconnected parts of the graphs ( = 0) and that the two 
points loop Nn (6, ^c) Nr, (Cc, 6) is replaced in (1.26) by "iV„ (C;,, Cc) iV„ (^c, 6) " 
which has no pole and is zero in rotationnally invariant systems. 

Consequently, since we have no poles at coinciding ^'s, the expansion in 
^ for the contribution of the p\ chains is the same whatever sector we choose 

Cai < ^ ••• ^ and the results must contain terms of the type ( |- 1 
with Jj > if p > 1. According to the chosen sector we use the convergent 
expansion 

1 1 ^ (i, 



Now, each chain separetely can be expanded in the choosen sector but if we find 
in the expansion some contributions ( ^ 1 with a negative or nul power Jj, 
we may ignore that contribution since we know that it must be cancelled when 
we sum over the p! chains. 

For instance, we consider three consecutive points of the chain and the cor- 
responding situation 

(Ci-1,6) (^.Ci+l) , ii = Inf{^i-,,^i,^i+l) (A.3) 
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according to (A. 2) we obtain the contribution 

This term provides for the variable a contribution 

-^Ji = -r-s-l<0 (A.5) 



so that this chain should be ignored in the choosen sector. In such a sector, 
only the chains with no other relative minimum that the end-points of the chain 
should be considered. 

A first important consequence of this remark is that the chains which form a 
loop gives no contribution; this can be verified for instance in (A.l) with a = d 

in the case of a three points loop. A second consequence is that if we compute 

-An{vt,^b,)l[^ni^b,,^b,+,) i^n(e6„%*) (A.6) 

i=l 

in the sector < ^ •■• ^ Cap then , is either or and along the 
chain we have 

61 < 62 < - < 6, = Cap > > ^6,+^ > - > 6, l<q<P (A.7) 
Prom (A. 2), we obtain for (A.6) 

-j^ / 1 \ / 1 ^ ri-i—ri + l 



n 



(A.8 



and by identification, we have 



which implies 



Jhi = i — n — r\ — \ (A. 9. a) 

Jb, = Ti-i -n 2<i<q-l (A.9.b) 

Jb, = rg-i + rg + 1 (A.9.C) 

A =Ti- n-i q + l<i<p-l (A.9.d) 

Jbp=n-j- rp-i (A.9.e) 

p 

J = J2Jb.=i-j (A.IO) 
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Clearly, the variables i and are all determined as function of n, j and Jt,. ; 
the constraints i > n + 1 and > together with Jf,- > implies only two 
constraints on j 

j<n-J+ (A.ll.a) 
j>n-J+-\ + l (A.ll.b) 

with 

p 

J+= J2 

The original sum over j of the function describing the kernels An (r?* , ) and 
(Cfcpi^l) is, for that chain, transformed into 

(A.13) 

where the operator is defined in (3.16). 

Finally, summing over all the p\ chains consists in summing over all subsets 

of {ai, ap_i}(including the empty set and the complete set) corresponding 
to the definition of J+. Taking into account the alternate sign in (— )^~^ the 
total sum is nothing but 

p n , 

11(1-^^.) E%^- (A-14) 
This result is the purpose of this appendix. 

B Appendix: the 't Hooft expansion 

We define the symbol 

n*-' = {n+l){n + 2)...{n + J) (B.l) 

which is expanded in powers of n as 

J 

n*^ = ^ (Tp ( J) n-^-P (B.2) 

p=0 



with 



ao ( J) = 1 aj ( J) = J! (B.3.a) 

(71 (J) = ■^(■^+^) C7p(J)=0 ifp>J (B.3.b) 
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For a given p, the coefficients ap (J) are polynomials in J; clearly we have 

(^P (J) = ^ k2...kp (B.4) 

l<fel<fe2<---<fcp< J 

and we observe the relation 

ap{J)^ap{J-l) + J ap-i{J-l) (B.5) 
Since ap (J) vanishes if J < p, we write 

(J) = (-^ - 1) - - P + 1) {J) (B.6.a) 

TTO (J) = 1 TTl ( J) = J + 1 (B.6.b) 

The relation 

J TTp (J) = ( J - p) TTp (J - 1) + 2pJ TTp-l (J - 1) (B.7) 

shows that Hp (J) is a monic polynomial of degree p in J which vanishes at 
J = -1 if p > O.We find 

7r2(J) = (J+l) 7r3(J) = J(J+1)' (B.8.a) 

7r4(J) = (J + l)(^j3 + j2-^J- (B.8.b) 

7r5(J) = J(J+lf (^j2-ij-2^ (B.8.C) 

We note that at large J 

p! \2 



^piJ)--y {-) (B-9) 



so that the so called "BMN" limit [n and J — > oo, ^ = cte) is 

n*-' ~ n-' e^^ (B.IO) 

From the generating functional 

C30 J 

(l + a;)-^-^ = ^(-)-^/'^^ (B.ll) 

and from 

^ (1 + xr-' = l^Al±^i (B.12) 
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we obtain by identification of the powers of x 

J+1 



3=0 



In section 4 we have to calculate 



j=JV-J 

From 

J 

{N - = iV (iV + 1) ... (iV + J) = ^ ap (J) TV-^-P+i 

(B.lS.a) 

{N-J- = TV (TV - 1) ... (AT _ J) = ^ ap (J) TV-'-f+i 

(B.lS.b) 

we obtain the topological t'Hooft expansion 

JV-l J 
j=N—J p odd 

Also, from (B.9), the BMN limit (large N, large J, = cte) is 

iV-l 



J V2A^ 



We find convenient to introduce the expression 



(B.lS.b) 

which has clearly a expansion 



„**(J+i)= ^ (-f S2p(J) n-^-2f+i (B.19) 

p=0 



27 



with 

So (J) = 1 
We have the relation 



S2 (J) = 



J(J+l)(J + 2) 
24 



J2 



^2p (J) = S2p ( J - 2) + — E2p_2 (J - 2) 

We write for p > 

y (J-2ff+2)(J-2ff+3)...J(J+l) 
^2p (^j (24f"^! 

^0 (J) = 1 ^2 ( J) = J + 2 

the relation (B.21) becomes 



(B.20) 
(B.21) 

(B.22.a) 
(B.22.b) 



J ( J + 1) (J) = ( J - 2p) ( J - 2p + 1) ^2p ( J - 2) + 6pj2 ^2^-2 (J - 2) 

(B.23) 

The polynomials 92p {J) are monic polynomials of degree p in J which vanish 
for J = -2 if p > 0. We have 



(J) = (J + 2) J + 



12 



26 ^ 48 



0,{J) = {J + 2)(J' + -J+- 

/ / . „N / 42 . 4184 ^ 576 

,(J) = (J + 2)(J^ + -J^ + ^J+- 



At large J 

so that the "BMN" limit (and the first correction) of n**(-^+^) is 



(B.24.a) 
(B.24.b) 
(B.24.C) 

(B.25) 



**(j+i) j+i 



e 24 



l + OI- 



72 

— = cie (B.26) 

n 



the exponential term being the first correction to the "BMN" limit. 
With these notations, we may compute directly the expansion of 



N-l 



fe=0 



r(fc + j + 1) 
r(fc + i) 



2P 



J- 



in 



jy**(7+l) 



where X^Li -^i = J'^ the "BMN" limit is trivially 



(B.27) 



r(fc + j+i) _ 2p A^-Q 



i=l 



fe=0 



r(fc + i) 



J 



2N 



JiJ_ 
N 



cte 



(B.28) 
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Finally, if we compare the coefficients of the term in the expansions 

(1 + x)-^-^ = (1 + x)-^-"-^ (1 + xf (B.29) 
we obtain the relation 

r(j + i) - ^ ^ ^ ^-^"^^ r(i + a + i) ^^-''^ 

if a is a non negative integer (with trivial analytic continuation to any a). As a 
consequence 

^ T{j + a + c+l)T{j + b + c+ 1) 

ho ^^^■ + ^) r(j + c + i) 



Inf{a,b) I ^„ 

'b 

(a + 6 + c-g+l) r(iV) 



g! C9 C« T{N + a + b + c-q+l) 



9=0 

(B.31) 
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